Interface steps in field effect devices 



m 

O 
O 
(N 

o 

Q 



o 
o 

u 

Oh' 



I 

O 

o 



o 

-I— > 



O 

o 



X 



Samuel WehrlQ and Christian Helm 
Theoretische Physik, ETH-Honggerberg, CH-8093 Zurich, Switzerland 
(Dated: February 2, 2008) 

The charge doped into a semiconductor in a field effect transistor (FET) is generally confined to 
the interface of the semiconductor. A planar step at the interface causes a potential drop due to 
the strong electric field of the FET, which in turn is screened by the doped carriers. We analyze 
the dipolar electronic structure of a single step in the Thomas-Fermi approximation and find that 
the transmission coefficient through the step is exponentially suppressed by the electric field and 
the induced carrier density as well as by the step height. In addition, the field enhancement at the 
step edge can facilitate the electric breakthrough of the insulating layer. We suggest that these two 
effects may lead to severe problems when engineering FET devices with very high doping. On the 
other hand steps can give rise to interesting physics in superconducting FETs by forming weak links 
and potentially creating atomic size Josephson junctions. 

PACS numbers: 73.40.Q, 74.78.-w, 85.30.Tv 



I. INTRODUCTION 

The field effect transistor is a widely used device, both 
for commercial products as well as in research on cor- 
related electron system (a recent overview is given in 
Ref. 0). Traditionally this technique is applied to semi- 
conductors, such as silicon or GaAs, and more recently 
to high temperature superconductors (see Ref. Ill2l and 
references therein) as well as organic materials^i^i^iSiL&. 
The field effect is intriguing because it allows in princi- 
ple to tune the charge density continuously by changing 
the applied voltage between gate and source/drain (see 
Fig. However, field effect experiments on the sys- 
tems mentioned above are difficult, because high electric 
fields, limited by the electric breakdown of the insulat- 
ing dielectric, are required in order to achieve substan- 
tial changes of the charge distribution. Nevertheless, new 
gate insulator materials, such as complex oxide dielectrics 
and ferroelectric oxides, in recent years allowed to push 
induced surface charge densities to promising valuesi. 
Thereby operating a FET at high electric fields leads to 
a strong confinement of the induced surface charge to the 
dielectric-semiconductor interface. Therefore, imperfec- 
tions of the interface become important. 

In the present paper the effect of interface steps in such 
devices is considered. Steps occur as imperfections when 
growing crystals*^, or can be created artificially. Natu- 
rally, one expects that steps act as barriers in the trans- 
port channel confined to the surface. Below we show 
under what circumstances this effect should be impor- 
tant, and that the transmission through the step depends 
exponentially on the electric field and the step height. 
Therefore, at high enough doping, a step cuts the trans- 
port channel and may present a severe problem. This 
may be particularly important in organic FETs where the 
step height is given be the molecular size and therefore is 
quite big. However, in superconducting FETs the step- 
barrier can act as a weak link and give rise to a Josephson 
junction where the critical current depends exponentially 
on the applied field. This might be an interesting way to 
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design a dissipationfree switch in a superconductor. 



We distinguish two types of field effect transistors. The 
term metal-insulator-semiconductor FET (MISFET) is 
used for a whole class of devices where the intrinsic car- 
rier density of the semiconductor is negligible [Ni k, 0) in 
the absence of the electric field. The electric field serves 
to induce a conducting space-charge layer on the surface 
(see Fig. nj. In contrast to this, in the superconducting 
FET (SUFET), the semiconductor is replaced by a su- 
perconductor or metal (above T^) which has a finite den- 
sity Ni. In this case the electric field alters the density 
which, for example, may change the superconducting T,.. 
In both cases the electric field induces a surface charge 
which screens the field. However, the extension zq of the 
surface charge perpendicular to the interface differs for 
MISFET and SUFET due to the diflFerence in Ni. As a 
surface step of height h is expected to be important for 
h ^ zq, we first give a rough estimate of zq for both cases. 



The problem of the charge profile in MISFETs was ex- 
tensively studied in Ref. II for the case of a continuous 
medium. Calculations for molecular crystals, which ex- 
plicitly take into account the discreteness of the lattice, 
were done recentlj*i2ii^. Below we follow Ref. ^3 and 
refer to this work for more details. In the case of a per- 
fectly flat and infinite interface, the wavefunction of the 
carriers can be separated into a plane-wave part paral- 
lel to the interface and a transverse part C(^)7 where the 
z-direction is perpendicular to the interface. The charge 
profile is then given by n(z) cx C(^)^- The transverse 
function is best calculated using a trial wavefunction 
for the lowest subband, such as the Fang-Howard trial 
wavefunction ({z) cx 2exp(— 6z/2), where b is the varia- 
tional parameter. We use this ansatz which is sufficient 
for our purpose, although more accurate trial functions 
existi^. Minimizing the energy with respect to b yields 
the average distance zq of the charge distribution from 
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the interface and the width w of the distribution: 
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In the equation above Ss is the dielectric constant of the 
semiconductor, its effective mass in the z-direction 
and no the induced surface carrier density. Note that we 
use a small n for surface density and a capital N for the 
volume density. Estimates of zo are given in Tab. U Pa- 
rameters for the Si-FET are standard and given in the 
literature. At the highest possible fields, zo = 12 A, 
which corresponds to a few unit cells. Pentacene was 
chosen as an example for an organic MISFET, and the 
parameters for the thin film transistor presented in Ref.0 
were taken. Only the effective mass is not accessi- 
ble by experiment and has to be estimated from band- 
structure calculations. Theoretical calculations yield a 
hopping tz = 0.47 meV between the Pentacene molecules 
in the z-directioni^. Assuming a tight-binding band, the 
effective mass at the minima of the band is given by 
nriz — / {2 a? tz), where a = 15.5 A is the layer spacing 
in the z-direction*. This yields — 0.34 mg. Equa- 
tion J^l gives then zq = 20 A which is of the same order 
as the layer spacing. Hence, molecular steps, as are ob- 
served when growing films, may be important in such 
devices*. Note that up to now, pentacene samples used 
in the experiments have many traps which may lead to 
corrections to the estimation of zo in Eq. (^). However, 
there are experimental attempts to use single crystals in 
order to reduce the number of traps^. 

In the case of SUFETs the superconductor has a finite 
density Ni of carriers which give raise to metallic screen- 
ing. Electric fields are screened within the Thomas- Fermi 
screening length Atf- Therefore zq « Atf for SUFETs. 
In cuprates (e.g. YBa2Cu307) Ni « 1-5 x lO^^ cm^^ 
and the Thomas-Fermi screening length is estimated to 
be 5-10 A (Ref. 2). This is larger than the width of the 
superconducting layer (3 A), but smaller then the unit 
cell (a = 15 A), which means that only the first supercon- 
ducting layer is affected by the electric field. Therefore, 
it is important that the transport in the SUFET occurs 
only in the first layer. Note that in a SUFET the drain 
and source electrodes usually contact the first layer (see 
Fig. However, the Josephson coupling among layers 
leads to a shortcut through the bulk. There are differ- 
ent possibilities to resolve this problem: (a) one uses a 
very thin superconductor which, ideally, is only one unit 
cell thick^^. (b) The distance between superconduct- 
ing layers is increased in order to suppress the Joseph- 
son coupling!—. (^Q^ T\ie interlayer Josephson coupling is 
suppressed e.g. by a parallel magnetic field (Fraunhofer 
like pattern) . In Ref. [l^ possibility (a) was realized and 
changes in Tc as well as an insulating phase were induced 
by the field effect. The ferroelectric Pb(Zra;Tii_3:)03 was 
used as a gate dielectric and a surface carrier density 
no = 7 X 10^^ cm~^ was achieved. The superconduc- 
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FIG. 1: Schematic picture of field effect transistor (FET). 



tor consisted of 1 to 2 unit cells of GdBa2Cu307_a;. 
Similarly, a single superconducting Cu02-layer has been 
created due to inhomogeneous oxygen doping (resulting 
in a distribution of Tc in the layers) on the surface of 
a Bi2Sr2CaCu2 08 single crystal in Ref. 0. Again, in 
this device the physics is confined to a small region at 
the interface and steps should play an important role. In 
particular, when the sample is superconducting, a step 
may induce a weak link. This effect can be of particular 
relevance in highly anisotropic layered superconductors 
such as the high-Tc material Bi2Sr2CaCu208 or inter- 
calated (LaSe)i.i4NbSe2, see Ref. 0, where the thick- 
ness/coherence length « 3 A of the superconducting 
(Cu02 or NbSe2)-layers is much smaller than their dis- 
tance (c? w 15 A), i.e. atomic-size steps can naturally 
constitute a weak link (case b) . 

In the following we solve first the problem of an inter- 
face step of a perfect conductor exposed to a perpendicu- 
lar electric field. This is a good approximation for h ^ zq 
and yields an induced dipole. Knowing the charge dipole, 
the potential barrier across the step is calculated using 
the Thomas-Fermi approximation. Finally, the transmis- 
sion coefficient for tunneling through the step is calcu- 
lated in the WKB approximation. 



II. PERFECT CONDUCTOR APPROXIMATION 

In the geometry with a step of height h along the y- 
direction (see Fig. O we model the interface by a two- 
dimensional grounded conductor with a potential differ- 
ence Vg to the gate in distance d. This approximation be- 
comes exact, if the thickness of the charged area is small 
compared to the step height, zq <ti h. Far from the step 
the electric field is homogenous, Eq = Vq/cJ « Vq/ (d+h), 
but it is distorted near the step, as a surface charge is 
induced there in order to compensate the potential dif- 
ference $0 = GVoh/d, which is realized for a flat inter- 
face with a homogeneous charge distribution. The typical 
equipotential curves near the step are shown in Fig.Q. 
As the system is invariant in the y-direction, we can use 
the conformal transformation (in the xz-plane) 



u = X + iz ■ 
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TABLE I: Comparison between diflerent FETs. and Sa 
are the dielectric constants of the dielectric and the semi- 
conductor respectively, is the effective mass in the z- 
direction which enters equation Q. A'i is the intrinsic car- 
rier density, Eh is the sample dependent breakdown field and 
no = i?i,ed/(47r e) is the maximum surface carrier density, zq 
is the distance of the charge distribution from the interface, 
which is calculated by Eq. Q for the MISFET and equal 
to the Thomas-Fermi screening length Atf in the case of 
the SUFET. h is the step height which, for Pentacene and 
GdBa2Cu307_2:, was taken to be the layer spacing in the z- 
direction. rii = Nih is the intrinsic surface carrier density 
in a single layer, kp = ^J2n{ni + no) is the Fermi- wave vec- 
tor. (T) is the average transmission coefficient calculated by 
equation "Ref.Q 'Ref.H =Ref.[ll ''RehEl "^See text, 
■''Ref.13 




FIG. 2; Geometry of an interface step. 



which maps the upper half complex plane onto the do- 
main above the step, the cuts being on the positive real 
axis. The solution for the potential is then given by 
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Far from the step, 'w{u) from equation ^ can be ex- 
panded in orders of h / y/ x'^ z^ , which yields the ap- 
proximate expression 

<^{x, z) = -Eo (^z-h ^^^^ , Vx^ + > , (5) 

where 9 is the polar angle in the xz-plane. In general, 
the surface charge of a perfect conductor is given by 




FIG. 3: Equipotential curves and electric field vectors near a 
step as in Fig.|5| 



n = edE±/{4:TTe), where E± is the electric field being 
perpendicular to the interface and Sd is the dielectric 
constant of the gate insulator. Using equation jS)) for 
the potential yields the asymptotic charge distribution 

{x > h) 
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with the doping no = £dEo/{4=TTe). The exact solution is 
given by the parametric expression 



n{x, z) = no 



1 — w{u) 



1 -I- w{u) 



(7) 



where w{u) is given by Eq. and u = x -\- iz \b, taken 
on the surface of the conductor. The exact result as well 
as the approximation equation © are shown in Fig 01 
Obviously, the electric field induces a dipole center at 
the step which falls off as x^^ far from the step. 

It can be shown by the exact solution 10 that there 



is a weak divergence 
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tively) of the surface charge at the upper corner of the 
step due to the sharp edge, which enhances the local 
tunneling rate through the dielectric and can therefore 
serve as a nucleation center for a possible breakthrough 
of the device. The singularity at the edge can be regu- 
larized: (a) at distances from the step smaller than the 
local Atf , where the quantum mechanical exclusion prin- 
ciple comes into play, (b) geometrically due to finite step 
curvature. In both cases the field enhancement is of the 
order {h/XcY^^ where Ac is the cutoff (Atf or finite cur- 
vature) which is at least of the order of the atomic scale 
(Ac « 1 — 2 A). Hence, ratios /i/Ac might be of the order 
10 — 20 and the field enhancement is expected to be a 
factor 2 — 3. This is not a big factor, but it should never- 
theless be important, because steps are line defects which 
have a big probability to hit a "weak spot" in the dielec- 
tric where an electric breakthrough can occur. Note also 
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FIG. 4: Surface charge on the perfect conductor (in units 
of no) along the step (distances in units of the step height 
h). Solid line: Exact solution with the divergences x~^^^ 
or (ft — z)~^^^ at (x,z) = (0,h). Dashed line: Approximate 
solution with the ~ 1/a; divergence. 



that kinks in the step would lead to a further enhance- 
ment of the electric field, (c) The singularity can also be 
regularized due to the finite bandwidth, where the locally 
induced charge exceeds the filling of the band. This oc- 
curs when the potential diff'erence $0 = eVoh/d exceeds 
the bandwidth. For the maximal gate fields ^0 = Ebh 
which is 2 eV for the Si-MISFET and 0.5 eV for the pen- 
tacene MISFET (see Tab. IJ). In the latter case this is 
roughly the bandwidth and a complete local charge de- 
pletion can be expectedi^. 



III. TRANSMISSION ACROSS THE STEP 

As discussed in the introduction, the wavefunctions 
of the electrons confined to a fiat and infinite interface 
can be separated into a parallel and a transverse part: 
4>kj{j^,z) = e^^'^Cji^) where r = {x,z) and k = {kx,kz) 
represent coordinate and wavevector parallel to the in- 
terface and j = 0, 1, 2 . . . labels the discrete and confined 
states in the transverse direction. The energies are given 



by 
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It can be shown that the eigenvalues ej of the tran- 

verse part scale as ej cx Eq^"^ cx Uq^^, where Eq and 
uq — SdEo I (47re) are the electric field and the surface car- 
rier density respectively (see Ref . Hi]) . If only the :/ = 
state is occupied (i.e. the lowest subband), then the sys- 
tem is two-dimensional and the Fermi energy is given by 
tp = no/g, where g = m||/(7r/i^) is the constant density 
of states. The condition for the occupation of only the 
lowest subband is < ei — eg which holds for densi- 
ties no below a certain threshold, in Si, this threshold is 



3 X 10^"^ cm~^ (using the triangular potential approxima- 
tion, see Ref. uM which is of the same order as the max- 
imally achievable surface densities (see Tab. IJ). Hence, 
it is not clear, whether only the lowest subband is oc- 
cupied. However, the presence of a step leads to charge 
dipole which suppresses the carrier density on the lower 
side of the step where the condition is then clearly satis- 
fied. As discussed below, it is exactly this region which 
is of interest when calculating the transmission coeffi- 
cient through the step. Therefore, we consider in the fol- 
lowing only the lowest subband and treat the system as 
two-dimensional. The non-uniform charge distribution, 
caused by the step, can be treated in the Thomas-Fermi 
approximation (TFA) which locally assumes a free (here 
2D) electron gas with density 



n(r) = gepir) = g[eF 



(9) 



where eF(r) is the spatially dependent Fermi Energy. 
At this point the result of the previous section enters 
by assuming that the surface carrier density n(r) = 
/ dz N{r, z) of the real system is still well approximated 
by the perfect conductor approximation (PCA). Note 
that in the PCA the induced charge is a pure surface 
charge with zero width and therefore differs from the 
real charge distribution which has a finite extension in 
the z-direction. However, the PCA becomes exact far 
from the step. Furthermore, knowing the solution of the 
PCA allows to calculate the local potential Vioc via re- 
lation In the following the approximate solution ® 
of the PCA is used which has the right asymptotics and 
which permits an analytic calculation of the transmission 
coefficient. Combining relations (0 and © yields 



V\oc{x) = = 

i: g X 
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The potential H10|l acts as a potential barrier across the 
step for a; > 0. In the WKB approximation the trans- 
mission coefficient T(k) of a mode k is given by 

T(k) = exp (^-2^ " ^{2nm/h^)mocix)-iu]dx^ , 

(11) 

where ek — Cko— eo— ky/{2nm) — fc^/(2m||) is the ki- 
netic energy along the a;-direction. Using the local poten- 
tial H10|) implies xo = noh/ijigt) = 2noh/k'^ and yields 
the analytic result 



T'(k) = exp — 27rno 



(12) 



Not surprisingly, the transmission coefficient depends ex- 
ponentially on the step height h and on the electric field 
which is proportional to the surface density no- In or- 
der to calculate the conductance through the step, the 
average (T) over the Fermi surface of the transmission 
coefficient enters, which for the parabolic band is 



{T) = J^^ ^T{kp cos ^) ^ f (^27Tno^^ , 



(13) 
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where the function / is given by 

2 /"^/^ , p^^ 

f{x) = - d^e--/-^^«W-^. (14) 

The approximate expression is vaUd for x 3> 1. In a 2D 
system with paraboHc dispersion the Fermi wave-vector 
is kp — ^27r(ni + no), where rii is the intrinsic carrier 
density. In the MISFET rii — which yields the simple 
expression (T) = f{kph). The conductance through the 
step is given by the Landauer Formula 

G = G„N{T), (15) 

where Go = 2e^/h w 12.9 kfl is the quantum conduc- 
tance. N = {2/tt) kpL is the number of channels for a 
sample of size L (along the j/-direction) taking into ac- 
count both spins and the condition —kp < ky < kp. 
The functional behavior of / implies that the conduc- 
tance through steps with height h > kp/{2TTnQ) (in the 
MISFET h > l/kp) are exponentially suppressed. 

Equation H13|l was applied to the FET examples dis- 
cussed previously and the results are summarized in 
TabQ] An arbitrary step height /i = 20 A, which might 
be due to an artificial step, was taken for the Si-MISFET, 
which yields (T) = 0.06 for the highest fields. In the Pen- 
tacene MISFET where l/kp 15 A which is almost the 
same as the interlayer spacing (a = 15.5 A) and which 
yields (T) = 0.2 for a molecular step {h = a) at the 
highest achievable dopings. 

Recently, there were unconfirmed claims that pen- 
tacene would become superconducting when doped to 
half fillingSS. In the light of our discussion such high dop- 
ings are unlikely, as the locally enhanced electric field at 
the step edges might lead to the electric breakthrough 
in the dielectric. Even if such electric fields could be 
reached in reality, the resulting device would be signifi- 
cantly limited by molecular steps, as for required density 
no = 2 X 10^"* cm~^ corresponding to l/kp = 3 A the 
transmission (T) = 2 x 10"'^ of the step would be small. 

In lower electric fields, for the SUFET below Tc, the 
charge dipole due to the step can form a Josephson junc- 
tion, which is of the SIS or SNS-type depending on the 
charge density in the nonsupcrconduting region. The 



critical current Ic through this junction can be estimated 
for a SNS tunnel junction by the Ambegaokar-Baratoff 
formula— 

/.i?./.G- = ^tanh(A,), (16) 

where A is the (s-wave) gap of the superconductor and 
G the conductivity from equation (|15|l . In this case the 
critical current Ic depends on the external electric field 
exponentially via the effective thickness of the insulating 
region, which enters the tunnel matrix element T. A sim- 
ilar exponential dependence Ic ^ exp^—l/lp) is found in 
field effect doped SNS-junctions with a large distance I of 
the superconducting leads, I ^ Ip — [hD / 2-^ ksT)^ 1/2) 
[D ~ n^/^ Diffusion constant in normal metal), while in 
the opposite limit / <C It the critical current depends al- 
gebraically on I, /c ~ 1/Z, see Ref. In both cases it 
is seen that the superconducing transport depends sensi- 
tively on the local charge density near the step, i.e. can 
be easily modified by an external electric field. 

Note that in a highly anisotropic layered superconduc- 
tors, such as high-Tc-materials an atomic size step can 
already form a weak link in the absence of the charge 
dipole effect described here, because the coherence length 
(f = 3 A) of the quasi 2D superconducting layers is 
smaller than the height h of the step. 

To conclude, we showed that a charge dipole is induced 
in the 2D electron gas near an interface step, which can 
lead to a (Josephson) weak link due to the local depletion 
of the charge density. This forms the limiting factor for 
transport through an ultrathin metallic or superconduct- 
ing layer and might be used for ultrasmall dissipationfree 
switches (SUFET). In addition to this, the field enhance- 
ment near the step in the insulating barrier can trigger 
the breakthrough of the dielectric and thereby limit the 
maximal doping level no, max- Both effects are hard to 
avoid and pose a fundamental challenge in term of atom- 
ically flat interfaces for any FET device with a quasi 2D 
(super)conducting charge density. 

We thank M. Sigrist, T.M. Rice and G. Blatter for 
useful discussions and the NRCC MaNEP of the Swiss 
NSF for financial support. 



* Electronic address: |swehrli@phys . et hz . cli| 
^ C.H. Ahn, J.-M. Triscone, J. Mannhart, Nature 424, 
(2003) 1015. 

^ J. Mannhart, Supercond. Sci. Technol. 9, (1996) 49. 

3 Y.S. Yang, S.H. Kim, J.I. Lee, H.Y. Chu, L.M. Do, H. Lee, 

J. Oh, T. Zyung, M.K. Ryu, M.S. Jang, Appl. Phys. Lett. 

80, (2002) 1595. 

Y.Y. Lin, D.J. Gundlach, S.F. Nelson, T.N. Jackson, IEEE 
Trans. Electron Devices 44, (1997) 1325. 
^ V.Y. Butko, X. Chi, D.V. Lang, A.R Ramirez, 
cond-mat/0305402 (2003); V.Y. Butko, X.Chi, 
A.P.Ramirez, . cond-mat/0307372, (2003). 



^ V. Podzorov, V.M. Pudalov, and M.E. Gershenson, 

Appl. Phys. Lett. 82, (2003) 1739; V. Podzorov, S. E. 

Sysoev, E. Loginova, V.M. Pudalov, M.E. Gershenson, 

cond-mat/0306192 (2003). 
^ R.W.I, de Boer, A.F. Morpurgo, T.M. Klapwijk, 

cond-mat/0307320 (2003) 
* J. Takeya, C. Goldmann, S. Haas, K.P. Pernstich, B. Ket- 

terer, B. Batlogg, cond-mat/0306206 (2003). 
^ M.J. Storcz, F.K. Wilhelm, Appl. Phys. Lett. 83 (2003) 

2387. 

° L.B. loffe, V.B. Geshkenbein, M.V. Feigel'man, A.L. 
Fauchere, G. Blatter, Nature 398, (1999) 679. 



6 



T. Ando, A.B. Fowler, F. Stern, Rev. Mod. Phys. 54, 
(1982) 437. 

S. Wehrli, D. Poilblanc, T.M. Rice, Eur. Phys. J. B 23, 
(2001) 345; S. Wehrli, E. Koch, M. Sigrist, to be published 
in PRE 

J. Sinova, J. Schliemann, A.S. Nunez, A.H. MacDonald, 
Phys. Rev. Lett. 87, (2001) 226802. 

R.H. Morf, N. d'Ambrumenil, S. Das Sarma, Phys. Rev. B 
66, (2002) 075408. 

J. Cornil, J. Ph. Calbert, J.L. Bredas, J. Am. Chem. Soc. 
123, (2001) 1250. 

C.H. Ahn, S. Gariglio, P. Paruch, T. Tybell, L. Antog- 
nazza, J.M. Triscone, Science 284, (1999) 1152 

F. Tafuri, J.R. Kirtley, P.G. Medaglia, P. Origiani, 

G. Balestrino, unpublished. 



M. Mosslo, R. Kleiner, R. Gatt, M. OneUion, P. Miiller, 

Physica C 341 (2000) 1571. 

P. Samucly, P. Szabo, J. Kacmarcik, A.G.M. Jansen, A. La- 
fond, A. Meerschaut, A. Briggs, Physica C 369 (2002) 61. 
M.R. Beasley (Chair), S. Datta, H. Kogelnik, H. Kroemer, 
D. Monroe: Report of the Investigation Commit- 
tee on the Possibility of Scientific Misconduct in the 
Work of Hendrik Schon and Coauthors, Sept. 2002, 
www. lucent . com/news_events/pdf /researchreview.pdf 
or publish.aps.org/reports/. 

V. Ambegaokar, A. Baratoff, Phys. Rev. Lett. 10, (1963) 
486 

A.F. Volkov, H. Takayanagi, Phys. Rev. B 53 (1996) 15162. 



